Kodama time: Geometrically preferred foliations of spherically symmetric spacetimes 
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In a general time-dependent (3+l)-dimensional spherically symmetric spacetime, the so-called 
Kodama vector is a naturally defined geometric quantity that is timelike outside the evolving horizon 
and so defines a preferred class of fiducial observers. However the Kodama vector does not by itself 
define any preferred notion of time. We first extract as much information as possible by invoking 
the "warped product" structure of spherically symmetric spacetime to study the Kodama vector, 
and the associated Kodama energy flux, in a coordinate independent manner. Using this formalism 
we construct a general class of conservation laws, generalizing Kodama's energy flux. 

We then demonstrate that a preferred time coordinate — which we shall call Kodama time — 
can be introduced by taking the additional step of applying the Clebsch decomposition theorem to 
the Kodama vector. We thus construct a geometrically preferred coordinate system for any time- 
dependent spherically symmetric spacetime, and explore its properties. We study the geometrically 
preferred flducial observers, and demonstrate that it is possible to deflne and calculate a generalized 
notion of surface gravity that is valid throughout the entire evolving spacetime. Furthermore, by 
building and suitably normalizing set of radial null geodesies, we can show that this generalized 
surface gravity passes several consistency tests and has a physically appropriate static limit. 

Keywords: Kodama vector, Kodama energy flux, spherical symmetry, Clebsch decomposi- 
tion, time-dependent metric, surface gravity, Hawking-Israel quasi-local mass, Hernandez-Misner 
quasi-local mass, Misner-Sharp quasi-local mass. Brown- York quasi-local mass. 



I. INTRODUCTION 

Black holes are an iconic part of Einstein's general rel- 
ativity. While we have a very detailed understanding 
of static and stationary black holes (the Schwarzschild, 
Reissner-Nordstrom, Kerr, and Kerr-Newman black 
holes), the situation with regard to evolving black holes, 
(evolving either due to accretion or Hawking radiation 
or both), is much more opaque. In particular, the rather 
limited number of currently known exact evolving solu- 
tions (Oppenheimer-Snyder collapse, the Vaidya solu- 
tion) makes it much more difhcult to fully describe an 
evolving black hole in any analytic detail. A fundamen- 
tal feature of the geometry of an evolving time-dependent 
spacetime is the lack of any (asymptotically timelike) 
Killing vector field, which seems to leave us without a 
preferred time coordinate with which to study the prob- 
lem. 

In 1980 Kodama made significant progress in this 
regard when he constructed a geometrically natural 
divergence-free preferred vector field that is guaranteed 
to exist in any time-dependent spherically symmetric 
spacetime This so-called "Kodama vector" defines 
a natural timelike direction in the region exterior to the 
black hole, and additionally induces an unexpected con- 
served current, but does not (in and of itself) define any 
naturally preferred time coordinate. By considering the 
"warped product" form of the spacetime metric for any 
spherically symmetric geometry we are able to investigate 
the Kodama vector (and the associated Kodama energy 
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fiux) in a coordinate independent manner. In particular 
we can use this formalism to easily generate a generalized 
Kodama fiux. 

By furthermore taking the extra step of invoking the 
Clebsch decomposition [2-4] on the (1+1) dimensional 
radial-temporal plane, we shall demonstrate that the Ko- 
dama vector field can indeed be used to construct a pre- 
ferred time coordinate, and more importantly a preferred 
coordinate system. The absence of any (asymptotically 
timelike) Killing vector in evolving spacetimes has made 
it difficult to achieve any consensus about best way to 
define such fundamental quantities as the surface grav- 
ity. Over the years, several different attempts have been 
made to extend the concept of surface gravity from static 
(and stationary) to time-dependent spacetimes. For in- 
stance, Hayward Q uses the Kodama vector itself as a 
substitute for the Killing vector, since it certainly pro- 
vides a preferred direction and it is parallel to the Killing 
vector in the static case (as well as at spatial infinity if one 
assumes the evolving spacetime is asymptotically fiat). 
Others (see 0] and references therein) have appealed 
to the freedom of normalization of the null geodesies to 
ensure their definitions reduce to known results in the 
static case Q. 

The layout of the current article is as follows: In section 
|TT]we briefiy summarize key properties of "warped prod- 
uct" spacetimes. In section Hill we present a quick review 
of the Kodama vector and Kodama's unexpected conser- 
vation law, and then significantly generalize Kodama's 
energy flux in section HVl Next, in section |Vl a Clebsch 
decomposition of the Kodama vector is made — in order 
to build a natural geometrically preferred coordinate sys- 
tem for any spherically symmetric time-dependent space- 
time. In sections IVII and I VIII we explore the Riemann 
and Einstein tensors in this geometrically preferred co- 
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ordinate system, being careful to connect the discussion 
back to the general "warped product" formalism of sec- 
tionHH Furthermore, in section lVIIII we review Kodama's 
conservation law in these preferred coordinates. In sec- 
tion IIXI we calculate the Brown- York quasi-local mass, 
and in section and |X] present our extended definition of 
surface gravity. Section IXII deals with the naturally in- 
duced notions of apparent and trapping horizon. Lastly, 
we add a brief discussion. 



II. WARPED PRODUCT SPACETIMES 

Any (possibly time-dependent) spherically symmetric 
metric can be written in the form 

ds^ = gab dx'^dx'' 

= %jix) dx'dx^ + rixf {dO"^ + sin^ d d^^} 

= %jix) dx'dx^ +r{xf "gc^fi dx^'dx^. (1) 

Here the two coordinates x^ run over the radial-temporal 
plane, while the two coordinates x" {9 and (f) run over 
the surfaces of spherical symmetry. The discussion can be 
generalized to {d+1) dimensions with d—\ dimensional 
spherical symmetry, but for now we are just working in 
(3-f 1) dimensions. (For higher-dimensional generaliza- 
tions in a Gauss-Bonnet context see I, Indepen- 
dent of the total dimensionality, gij is a (1-1-1) dimen- 
sional Lorentzian metric. 

Geometrically this is called a "warped product" man- 
ifold, with the radial-temporal plane being referred to 
as the "base space" , the surfaces of spherical symmetry 
being referred to as the "fibres" , and the function r{x) 
which depends only on the base space coordinates being 
referred to as the "warp factor" . It is a standard compu- 
tation to show that (up to the usual permutation sym- 
metries for the indices) the only non-zero components of 
the Riemann tensor (in any warped product spacetime) 
are 

Rijkl — ^Rijkl] (2) 

Riaji3 = - r {ViVj-r} ""gap; (3) 

RaPf^i^ = {"Ralltiiy " |Vrp {"gat^'^gflu " ''Saiy^'g/Jp) }• (4) 

An abstract computation along these lines can be found 
in O'Neill [TTj, page 210] while more explicit computa- 
tions can be found in |12h14| . Note that the covariant 
derivatives appearing above are covariant derivatives in 
the base space. But because r{x) depends only on the 
(1-1-1) dimensional base space coordinates, and because 
of the specific form of the warped product metric, these 
derivatives can be "bootstrapped" to covariant deriva- 
tives in the total warped product spacetime. 

In the specific situation we are interested in the base 
space is two dimensional, so in terms of the Ricci scalar 
of the radial-temporal plane we have the specific simpli- 
fication 



^R 



ijkl 



(5) 



Furthermore the fibre is a constant curvature sphere of 
radius unity, so 

"RaPiii' = Cgafj.'^giiL' - ''gau^gppi) ■ (6) 
Thus we now have on purely geometrical grounds 



Rijki — -Y i^k'^gji - %i''gjk) ; 



(7) 

RiajH = ~r {VjVjr} "gap; (8) 

Rcp^u^^r^ {l-|Vr|2} ("5aM"5/3- - • (9) 

It is also common in spherical symmetry to define 
the Hawking-Israel / Hernandez-Misner/ Misner-Sharp 
quasi-local mass O, [l^ by 



2™ 

1 = Vr ■ 

r 



(10) 



where both m{x'^) and r(x') are scalar functions on the 
base space. We now have 



Rijki - -j- {""gik^gji ~ %i"g3k) ; 



(11) 

R^c.JP = - r {V.Vj-r} "g^p; (12) 
Rapfiu ^ 2mr Cgafj.'^gpi, - '"gai^^gpii) ■ (13) 
It is often useful to go to an orthonormal basis, in which 



case 



^R / 

^ijici = 'Y [^ik ^ji ^ hi ^]k J ' 



R"- - '•■a 



(14) 
(15) 



R ■ = — 



(S&i, (5^- - 6^^ . (16) 
For the Ricci tensor we have 



Ria = 0; 

f 2m 2 1 F 



and for the Ricci scalar 



i? = ^i? - 4- 



. V r Am 



r.3 ' 



(17) 
(18) 

(19) 
(20) 



where the Laplacians above are in the (1-f-l) dimensional 
sense on the base space. The Einstein tensor takes the 
form 

2{V,Vjr} (Tsj'^r 2to 1 



Gia = 0: 



G. 



a/3 



%v (21) 
(22) 

(23) 



These are all purely geometrical statements — while one 
has chosen coordinates cc" = (x'; x°') to make the warped 
product structure manifest, these results are completely 
independent of one's choice of coordinates on the 
radial-temporal plane (the base space), and for that mat- 
ter are completely independent of one's choice of coordi- 
nates on the spherically symmetric fibres. 
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III. THE KODAMA MIRACLE 

It is well known that in a time-dependent spacetime, 
there is no (asymptotically timelike) Killing vector to de- 
fine a preferred time coordinate. The calculation of im- 
portant quantities such as the four-acceleration and the 
surface gravity become much more ambiguous. Addition- 
ally, there is no general consensus on the "best" form of 
the metric, nor on the "best" choice of the coordinate 
system. 

An interesting insight on this problem is given by Ko- 
dama [H, who proved the existence of a divergence- free 
vector field for any time-dependent spherically symmet- 
ric metric. The Kodama vector, fc", lies in the (1+1) di- 
mensional radial-temporal plane, so that k"^ = (fc'; 0,0). 
More precisely 



k" = 



(24) 



where the tensor is the (1-1-1) dimensional Levi-Civita 
tensor in the radial-temporal plane, denoted e^, canon- 
ically embedded into (3-f 1) dimensions according to the 
prescription 



ij 


o" 


. 


0_ 



(25) 



It is straightforward to check that Var — 0. Fur- 
thermore, if we define a positive semi-definite norm by 
||fc||2 = \g{k,k)\ = \g-\k\k% then ||fc|| = ||Vr||. (We 
shall use the superscripted symbol t> to denote the pro- 
cess of turning a vector into a covector by "lowering the 
index" , and use the superscripted symbol jl to denote the 
inverse process of turning a covector into a vector by 
"raising the index".) By appropriate choice of orienta- 
tion on the radial-temporal plane one can choose /c" to be 
(asymptotically) future pointing. It can also be defined 
by the more abstract statement 



k = (*2 dry 



(26) 



where by this one means "calculate the one-form dr, ap- 
ply the (1+1) dimensional radial-temporal Hodge star 
operation, and use the metric to convert the resulting 
one-form to a contravariant vector" . At this point it is 
necessary to emphasise, as originally pointed out by Ko- 
dama himself, that the Kodama vector does not in gen- 
eral reduce to the Killing vector in a static spacetime; all 
that one can say in general is that in static spacetimes it 
is parallel to the Killing vector. In regions where the Ko- 
dama vector is timelike (and we shall [informally at this 
stage] refer to this as the black hole exterior region, i.e., 
the domain of outer communication) the Kodama vector 
defines a preferred class of fiducial observers (FIDOs) [lH| 
specified by the unit timelike 4- vector 



V 



(27) 



Although the Kodama vector provides a preferred "time 
direction", and simplifies the evolution equations of a 



dynamical spherically symmetric system [16|, |l7[ , it does 
not at this stage define a preferred "time coordinate". 
We shall subsequently use the Kodama vector plus the 
the Clebsch decomposition theorem to construct a geo- 
metrically natural preferred time coordinate. 

To prove that the Kodama vector is divergence free the 
best strategy (with hindsight) is to consider the quantity 



--da 






o" 


_ 


0. 



;0,0 



^(-V..V;0,0 
0. 



(28) 



Note that the last covariant derivative is a base space co- 
variant derivative which vanishes since the (1+1) dimen- 
sional Levi-Civiata tensor is covariantly constant with re- 
spect to the (1+1) dimensional covariant derivative. But 
this (3+1) dimensional result, Va(e'j''/r^) — can easily 
be rearranged to give 



which now implies 



Vafc'^ = 0, 



(29) 



(30) 



so the Kodama vector itself is conserved. In addition, 
Kodama also proved that in time-dependent spherically 
symmetric spacetimes there is another (somewhat un- 
expected) conserved current. In terms of the Einstein 
tensor and the Kodama vector we have J"' = k^, 
with: 



VaJ"=V„(G"'' kb)=0 



(31) 



This is a purely geometrical statement, ultimately due to 
the warped product form of the metric — it is not related 
to the Bianchi identities. Let us specifically compute 

.r = G"'' h - (G'^kf, 0, 0). (32) 

Now working on the radial-temporal base space we have 

V, V,r f 2V2r 2m1 , ,, , , 
:r> k'. (33) 



Gij fc"* 



But using the fact that in (1+1) dimensions 
we have 



(34) 



V''r 



1 
2 



= +vV h 



Vfcr e±y} V'^r 



(35) 
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Combining these results 

G^jk^ = +^e±.,Vn|Vr|2} - ^ k\ (36) 

But in view of the definition of the Hawking-Israel/ 
Hernandez-Misner/ Misner-Sharp quasi-local mass [l^, 
il9| we then have 

1 2m 

G,, F = --e^^^yi{2m/r} - — k\ (37) 

That is 

G„- k^ = -^e^y V^m. (38) 

This computation has been performed using the (1+1) 
dimensional covariant derivative in the radial-temporal 
base space, but at this stage we can safely use the sym- 
metries of the situation to lift this equality to the full 
spacetime 

r = G"''fcb = ef Vfcm. (39) 

This is a purely geometrical statement — fundamentally 
connected with the warped product nature of the space- 
time — that underlies the unexpected conservation of 
the Kodama current. In view of the fact that we have 
already proven Va(e°''/r^) = we finally see 

VaJ" = ef Va^bm = 0. (40) 

Thus conservation of the Kodama flux is is a subtle result 
deeply connected with the warped product nature of the 
spacetime. We have presented this derivation in some 
detail because it is now possible to rapidly generalize the 
result in a significant manner. 

IV. GENERALIZED KODAMA FLUX 

Consider an arbitrary function $(TO,r) of the two 
quantities m(x^) and r(x*). Now construct the current 

Ji = {a„$(m,r) G^" -2dr^m,r) g'^"} ka. (41) 

This current is conserved in any spherically symmetric 
spacetime. To prove this note that by the definition of 
the Kodama vector and the geometrical identity proved 
above we have 

Jl = -2-^ {9m«'(TO, r) Vbrn + 5r$(m, r)V&r} 

= -2^Vb$(TO,r). (42) 

Conservation of this 4-vector is then obvious from the 
last expression. Note in particular that by the above 
argument any flux of the form 

^f2 = {/iMG"'' + /2(r)ff"''}A:a (43) 



[for arbitrary fi{m) and /2('')] will automatically be con- 
served. 

Formally, there is an even more general conserved cur- 
rent one can write down: For any arbitrary scalar func- 
tion ^'(a;*) defined on the radial-temporal base space the 
quantity 

= ^ Vh^, (44) 

is conserved. Though this result is more general, it 
is somewhat less geometrical, and does not have the 
same flavour as the above. If (and only if) the func- 
tions r{x^) and m(a;') are functionally independent (so 
that one can use m and r as coordinates on the radial- 
temporal base space) then these two notions ( J$ and Jip) 
can be made to coincide. In particular, ^ — >■ —2m{r,t) 
gives us Kodama's conserved flux J", while — >■ is 
just the statement that the Kodama vector itself is con- 
served, Va^" = 0. (For related comments in a higher- 
dimensional Gauss-Bonnet context see [ol. [lOj.) 

V. KODAMA TIME 

The Kodama vector has been used before in several as- 
pects of the time-dependent gravitational collapse prob- 
lem. However it has not been used to obtain a preferred 
time coordinate, nor a preferred coordinate system for 
the metric of a dynamic spacetime. Fortunately, in (1-f 1) 
dimensions it is possible to use the lesser known but clas- 
sic Clebsch decomposition theorem, a result complemen- 
tary to the more usual Helmholtz decomposition theo- 
rem, (see for instance 0-|3) to assert that there are two 
unique scalars a and j3 such that the Kodama covector 
k^ takes the form 

k^ = a dp. (45) 

Now in the "normal" exterior region where dr is space- 
like, (i.e., in the domain of outer communication), the 
Kodama vector and covector are both timelike, so in this 
region the one-form d/3 is guaranteed to be timelike. This 
very strongly suggests that /3 should be adopted as a pre- 
ferred "time coordinate". In fact, relabeling the scalar 
/3 as a coordinate (3 t, and relabeling the integrating 
factor as a — >■ i^, we write 

k^ dt. (46) 

This naturally induces a geometrically preferred time co- 
ordinate t — which we shall refer to as the Kodama time. 
Using this time coordinate is at least as natural as us- 
ing r for the radial coordinate. (This key step, though 
mathematically elementary, goes well beyond anything in 
Kodama's original paper [l|, or the various papers that 
have subsequently sought to use Kodama's formalism.) 
As we shall soon see, this choice of Kodama time coor- 
dinate is the unique choice that makes integral curves of 
the vector dt coincide with integral curves of the Kodama 
vector. That is 

k(xdt. (47) 
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Ultimately, adopting these coordinates (no matter how 
natural they appear) is of course a choice, and will be 
"justified" only insofar as they turn out to be useful. 

Adopting these (t, r) coordinates as preferred coordi- 
nates on the radial-temporal plane, and without any loss 
of generality, the metric can be written as 

ds^ = gtt{r, t) At^ + 2gtr{r, t)dr dt + g„{r, t) dr^ 

W {d6'2+sin2 6id<?!)2}.(48) 

However, since the Kodama vector k is orthogonal to dr, 
then also dt is orthogonal to dr, and so the cross term in 
the metric is zero. Thus in these preferred coordinates 
the metric is diagonal, 

ds^ = gtt(r,t)dt^ ^ grr(r,t)dr'^ ^r^ {dB^ ^ sin^ d d(i?) . 

(49) 

This relatively long argument has ultimately led us back 
to one of the simplest, and arguably most obvious, forms 
of the metric — a simple diagonal metric in Schwarzschild 
curvature coordinates. (Of course now we can argue that 
we have a geometrically natural reason for adopting this 
particular set of coordinates.) 

When using the Schwarzschild radial coordinate r it is 
natural to choose the parameterization 



grrir,t) = 1 - 



2r7!,(r, t) 



(50) 



for the radial-radial part of the metric tensor. Doing so 
will automatically give the quantity m{r,t) a natural in- 
terpretation in terms of the Hawking-Israel / Hernandez- 
Misner/ Misner-Sharp quasi- local mass [la, [l^. Since 
the radial-temporal plane by definition has Lorentzian 
signature, this choice then guarantees that it is possible 
to write the temporal-temporal component of the metric 
tensor in the form 



gtt{r,t) 



-2*(r,t) 



1 



2m(r,t) 



(51) 



We finally have the (quite standard) result 



= _e-2*('-^*) 1 _ 



2m(r, t) 



dt^ 



dr2 



r J 1 — 2m{r,t)jr 

-Kr^ {d6|2 +sin2 6'd(/)2} , (52) 

where in addition we know 

V^r = (dr), = (0,1; 0,0), (53) 

and the equivalent contravariant result 



(d.)'^^(0,l-^^;0,0 



(54) 

Furthermore, the components of the Kodama vector and 
covector in these coordinates are 

r = e*('''*)(l,0; 0,0); 

fc. = -e-*(-*)(^l-^^)(l,0;0,0). (55) 



As previously mentioned, the squared norm of the Ko- 
dama vector is equal to that of Vr: 



iVrl 



2m(r, i) 



(56) 



In these coordinates it is useful to define the time trans- 
lation vector, T, which is not a Killing vector unless the 
geometry happens to be static, as 



T^du 

-2*(r,t) , ^ 



r'^ = (l,0; 0,0) 
2m(r,i) 



(1,0; 0,0). (57) 



The squared norm of T is equal to the absolute value 
of the temporal-temporal component of the metric, \gu\i 
and proportional to the squared norm of the Kodama 
vector: 



That is, 

-2*(r,t) _ II'^IP . 



2TO(r, t) 



e-2*(-.*) ||fc||2. 

(58) 



imi 



(59) 



In the static situation it is the time translation vector 
T that will reduce to the (asymptotic) timelike Killing 
vector: T ^ K ^ k. Because the normalizations of T 
and k differ, then if the Kodama vector is used simply 
as a substitute for the Killing vector when attempting to 
calculate quantities such as the surface gravity [5] , one is 
likely to encounter normalization issues when taking the 
static limit. 

Specifically, to obtain a finite value of the four- 
acceleration near a possible horizon, as measured by an 
observer at infinity, it is necessary to multiply by a suit- 
ably defined normalizing factor [lO, [2l|. In the static 
case this normalizing factor is just \gu\^ and coincides 
with the squared norm of the Killing vector. However, 
in the time dependent case, not only does the geometry 
not possess a Killing vector, but also the squared norm 
of the Kodama vector does not coincide with \gtt\- This 
leaves us with a somewhat ambiguous situation with re- 
spect to the normalizing factor and the surface gravity of 
a time-dependent metric tensor, and means that we will 
have to exercise some care in defining the surface gravity 
of a time-dependent geometry. 



VI. RIEMANN TENSOR 

It is now a standard exercise to calculate the various 
components of the Riemann tensor (for instance, by using 
Maple). We note that the Riemann tensor is considerably 
less fearsome than one might suppose. Only one compo- 
nent is in any sense "difiicult" . Temporarily suppressing 
the (r, t) arguments for conciseness, and working in an 
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orthonormal basis we have: 



2m 



1 - 



Ot 



r 

2™ 



2m' 



[-$" + ($')'] +3(m/r)'$' 



(60) 



In view of the warped product formalism, we know that 
this rather messy quantity has a drect and simple physi- 
cal/ mathematical interpretation: As may be verified by 
direct computation it is simply ^i?/2, one half the Ricci 
scalar of the (1+1) dimensional radial-temporal plane. 
The remaining components are much simpler: 



R 



't(j>t<t> 



2m\ $' , , 

1 -; (61) 

r J r 



^t9f9 ^iifi r^{i-2m/ry 



(62) 



(63) 



These three quantities are easily seen to be proportional 
to ViVjT. (In fact they equal — {VjVjr}/r.) Finally 



R 



2m 



(64) 



Note that the particularly simple formula for R§^§^ un- 
derlies the identification of m{r,t) as the Hernandez- 
Misner quasi-local mass [l8l |. 



VII. EINSTEIN TENSOR 

For the Einstein tensor, working in an orthonormal 
basis, the single most important result is 



^tt — 



2m'{r,t) 



(65) 



Here the primes denote differentiation with respect to r, 
and the dots with respect to t. This result is utterly stan- 
dard, with the only novelty being that this formula for 
continues to hold in the time dependent case (subject 
of course to the coordinate choices made above) . This re- 
sult for is intimately related to the physical interpre- 
tation of m(r, t) as the Hawking-Israel quasi-local mass. 
A second important result is more subtle: 



^t> — 



2 TO(r, i) e*('''*) 



2m(r,f) 



(66) 



We shall soon see that this formula for G^- is central 
to the coordinate-based verification of Kodama's unex- 
pected conservation law, and that it is intimately related 
to the Brown- York quasi-local mass [22i] . 



For completeness we indicate 



2m'{r,t) 2<l>-(r,t) (l - ^) 



(67) 



This result is quite standard, (see for instance equation 
(2.65) of [1^), with the only novelty being that this for- 
mula for Gff continues to hold in the time dependent 
case. This now implies the useful result 



2$ 



G;; 



'ir,t) (l 



2m(r,t) 



(68) 



Finally, now suppressing the (r, t) arguments for con- 
ciseness, we have 



Ga 



Gi 



1 - 
3$' 



2m 
r 
m 

o2* 



1 



(69) 



m 
r 



m 
r 



1 



TO 

r 



4(to) 



<I> TO 



fl-2i») 



The first two lines here are again quite standard, and ap- 
pear also in static situations. (See for instance equation 
(2.66) of [1^). All the time derivatives have been iso- 
lated in the third line. With a little more work this can 
be somewhat tidied up as follows 



Gaa — G77 



TO 

r 



Tyjl — 2m/r 



Ot 



r 1 



o \ 3/2 



1 



(70) 



This rather complicated expression can be verified to 
equal — ^i?/2 + V'^r/r, the result we expect based on 
the warped product formalism. We note that the time 
derivative contributions to the Einstein tensor are quite 
isolated, and in this geometrically preferred coordinate 
system occur only in the Ggg = G^^ and G^~ compo- 
nents. This ultimately one of the key reasons we will 
find the Kodama time to be so useful. 



VIII. CO-ORDINATE BASED VERSION OF 
KODAMA'S CONSERVATION LAW 

Kodama's conservation law can now be studied in more 
explicit coordinate-based detail. First, based only on 
spherical symmetry and the definition of the Kodama 
vector, the unexpected conserved current J° takes the 
form 



J'^ = 5°'' Gbc k'= = {-fc^A:'' + (d^^)'^(d^^)''} Gbc y 



(71) 
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whence, since by construction ||fc|| = ||dr||, we see 

r = -Giik'^ + GiA'^rr- (72) 

But we have already exphcitly calculated the quantities 
Gii, Gif, k" and (dr)^ We obtain 



(73) 



This vector is is now obviously conserved since the 4- 
divergence is simply 



2 

m'(r, t) + m! [r, t)] 



g-<I'(r,t) y,2 

= 0. 



(74) 



Equivalently we note that from this coordinate-based cal- 
culation we explicitly recover 

-2 ft^. (75) 



This relation is somewhat miraculous in the present co- 
ordinate based calculation, and as we have seen has a 
deeper justification in terms of the warped product form 
of the spacetime geometry. 



IX. BROWN-YORK QUASI-LOCAL MASS 

With the coordinate system developed above, the no- 
tion of quasi-local internal energy arises naturally as the 
Brown- York quasi- local mass [22|. To prove this, first let 
us take some imaginary spherical surface surface r = rg, 
and hold tq fixed in time. Then the total energy inside 
this spherical surface depends on the net flux across the 
surface. To calculate the net flux we need the Gtr com- 
ponent of the Einstein tensor, in an orthonormal basis. 
That is 



Gif{r,t) 



2m(r, t) 



(76) 



r^(l — 2m(r, t)/r) 

whence, via the Einstein equations Stt Gab = Tat, we have 
the flux density 



m(r, t) 



(77) 

4 7r r2(l - 2m(r,i)/r) ^ ' 

Now the total net flux across the imaginary surface at 
r = To, in an amount of proper time r, is 

(net flux)j"'"' ^— I f{ro,t)x (area) x dr 



X dt. 



= / f{ro,t) X (An 4) 



2m{rQ,t) 



(78) 



That is: 

(net flux' 



ttinal 
tinitial 



■m{ro,t) 




ttinal 



(79) 



Then, if initially there is no mass inside r = rg, at time 
tfinai we have 



(net flux) 



tfinal 
tinitial 



'^o 1 



2 m(ro, tfinal) 



(80) 



In this situation, the only meaningful deflnition of inter- 
nal energy att = tinitial is to set U (ro, tinitial) = 0. Hence, 
at any subsequent time t the internal energy U{ro,t) is 
equal to the net incoming flux and so it makes sense to 
define 



2m(ro,t) 
ro 



(81) 



This internal energy is just the Brown- York quasi-local 
mass for the spacetime geometry with metric ([5^ . see 
[23 [. We can rearrange this (as pointed out in J2]) to 
yield 



m(ro,t) = U{rQ,t) - 



2ro 



(82) 



Here m(ro,t) retains its interpretation as the Hawking- 
Israel (and Hernandez-Misner/ Misner-Sharp [l^ [l^) 
quasi-local mass. The difference between the two notions 
of energy is just the self interacting Newtonian gravita- 
tional potential of a massive shell of radius tq. Both 
energies coincide at spatial infinity with the ADM mass. 



X. SURFACE GRAVITY 

Several attempts at calculating the surface gravity for 
a time-dependent metric have been made using the Ko- 
dama vector instead of the Killing vector , with results 
qualitatively similar to those in the static case; even to 
the extent of deriving some form of the first law of (black 
hole) thermodynamics. 



X.l. Surface gravity from fiducial observers 

The most intuitive way of calculating the surface grav- 
ity is by working in the exterior region and considering 
the four-velocity V parallel to the Kodama vector. Cal- 
culate the four-acceleration A = S/yV. Then explicitly 
computing the magnitude of this four-acceleration we see 



1 



2m(r,t) 



m{r,t) m'{r,t) 
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2m(r,t) ,, 
1 ^^$'(r,t). 



(83) 



(Note that near spatial infinity we have the sensible 
Newtonian result a — >■ m/r^.) The surface gravity can 
be defined as the acceleration of an observer near the 
evolving horizon, which we implicitly define by rnit) = 
2m(rH(t),t), as measured by another observer at infin- 
ity. Thus, at this point it is necessary to multiply by a 
normalizing factor, often referred to as a redshift factor. 
In the (asymptotically fiat) static case there is no doubt 
that the appropriate normalizing factor is 



\K\ 



|5,,|V2 ^ e-*W^i_2m(r)/r (84) 



and that the appropriate object to consider is the near 
horizon limit of 

/tstatic^ll^ll \\K\\^\\VkV\\ 



m{r) m'{r) 



2m(r) 



1 - 



$'(r) 



(85) 



The location of the horizon is in the static case implicitly 
defined by rn = 2m(rif ), and this is now a true Killing 
horizon and also an event horizon, at which we have the 
standard result M 



Kstaticlif = e ^ 



1 - 2m'{rH) 
2rH 



(86) 



When bootstrapping to the dynamic case a plausible 
generalization (which we shall subsequently buttress by 
also considering the radial null geodesies) is to replace 
1 1 if II — >■ ||T||, which at least has the virtue of maintain- 
ing the correct static limit. Under this hypothesis the 
appropriate object to consider is the near horizon limit 
of 

^v = P||||r|| = ||Vrl^|| 



^ g-*('-.t) 



rn(r, t) m'(r, i) 



1 - 



2m(r,t) 



$'(r,t) 



(87) 



which on the evolving horizon has the limit 



2rH{t) 



However, this proposed definition presents us with a po- 
tential ambiguity — what is the physically most appro- 
priate choice for the normalizing factor? The choice of 
I |T| I as normalizing factor as advocated above is certainly 
plausible, and has the correct static limit. Furthermore 
it is intimately related to the Kodama time introduced in 
this article, rather than the Kodama vector k. Neverthe- 
less, it is useful to see if we can come up with other plau- 
sible candidates for surface gravity in an evolving space- 
time, and see whether they agree with (or are closely 
related to) the above proposal, and whether they possess 
the correct static limit. 



X.2. Surface gravity from radial null geodesies 

There are other, possibly less ambiguous, ways to use- 
fully define the sought after surface gravity. More specif- 
ically, we can parameterize the strength of the gravita- 
tional field throughout the entire spacetime geometry by 
using the inaffinity properties of the radial null geodesies. 
Consider (temporarily) the following null vectors: 



(89) 



In the exterior region (where k is timelike) these null 
vectors are both outward-pointing, ^+ is future-pointing, 
and £~ is past-pointing. (Note that is then inward 
pointing; these specific conventions have been chosen 
to simplify the computations below as far as possible.) 
These are arguably the simplest radial null vectors one 
could write down using only the Kodama vector. It is 
easy to check that £°l = |||A:|p. Since we are working 
with spherical symmetry, both radial null vectors must 
satisfy the geodesic equation (in its non-affine parame- 
terized form): 



Vi^£± = K£± i± 



(90) 



where ki± are scalars defined everywhere throughout the 
spacetime. By contracting these equations with ^J, we 
can explicitly compute k^^ , to yield: 



m{r,t) m'{r,t) 1 



2 TO(r, t) 



(91) 

(Note that near spatial infinity we have kg^ — > m/r^.) At 
the evolving horizon this would reduce to the tentative 
definition 



knit) = 



l-2m'{rH{t),t) 
2rHit) 



(92) 



Unfortunately this does not reduce to the known re- 
sult in the static case — there is a missing factor of 
g-HrHit},t) _^ e-*(''"). This makes the above defini- 
tion not suitable for calculating the Hawking tempera- 
ture. (We emphasize that in static situations the stan- 
dard Wick-rotation trick of going to Euclidean signature, 
demanding the absence of any conical singularity at rn, 
and interpreting the Hawking temperature in terms of 
periodicity in imaginary time, uniquely enforces equa- 
tion ([55)) as the onlyphysically acceptable candidate for 
the surface gravity Q.) The source of the difficulty is, 
since e^*^''*^ = ll'7'||/||fc||) ultimately due to the fact that 
||T|| ^ ||fc|| in general. 

These considerations do suggest an improved strategy: 
Since we have seen how to use the Clebsch decomposi- 
tion to deduce the natural existence of a Kodama time, 
in addition to a Kodama vector, then it would seem ap- 
propriate to use the Kodama time as the natural (non- 
affine) parameter for these radial null curves. (That is, 
we now parameterize the null curves by Kodama time. 
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rather than the usual KiUing time used in the static case.) 
This is tantamount to choosing 



1 



TTfcF 



-*(r,t) ~n±. 



(93) 



This time, the inner product \s, l°L ^ \\\T\^ ■ These 
"Kodama time normahzed" radial null vectors are again 
tangent to the radial null geodesies and so satisfy 



A brief calculation yields 



whence 



Ki, = e 



(94) 
(95) 



1 



2m{r,t) 



$'(r) i Ti$(r,t). (96) 



At the evolving horizon, reduces to 

This is not quite equal to Ky^j^ — though it does share 
with Kv\H the desirable property of having the correct 
static limit. An improved proposal is to average over 
past and future pointing null geodesies and take 



(97) 



Knull = - + 



(98) 



Then KnuU = I'^v- That is: If one takes future-pointing 
and past-pointing outward null geodesies, normalized to 
Kodama time, and averages the resulting inaffinity pa- 
rameters, then one obtains the same Ky that we tenta- 
tively identified based on the 4-acceleration of the FIDOs 
that follow the Kodama flow. 

In short: By using Kodama time in addition to the 
Kodama vector we have now developed a geometrically 
preferred notion of surface gravity for spherically sym- 
metric evolving spacetimes that can meaningfully be ex- 
tended throughout the entire spacetime, and in addition 
exhibits a good static limit. 



XI. THE EVOLVING HORIZON 

With the calculations presented so far, it is not possi- 
ble to conclude too much about the evolving horizon at 
rnit) — 2 m{t, ru)- To relate this to a trapping horizon, 
in the Hay ward sense @| i it is necessary to compute 



the expansions, 9^^ and Qi^ , of the radial null vectors. 
Let us use the following definitions for the expansion 7] 



3t^=Vat±-tii^\ (99) 



A brief computation yields 



If 2m{r,t) 



(100) 



In particular both 9^;^ and 9i^ , change sign at the evolv- 
ing horizon ruit). This is sufficient to guarantee that the 
evolving horizon at r/f(<) = 2m{t,rH{t)) is an apparent 
horizon. However rnit) is not a trapping horizon in the 
Hay ward sense Q, nor an "Ashtekar horizon" [l^, since 
in addition 9^ and 9i_ also both vanish on the evolving 
horizon. This is not critical for our purposes, since ulti- 
mately a trapping horizon is not needed to have Hawking 
radiation |25l - [27| . 

For completeness, we have also computed £°L'Va9jj^ and 
£"LWa9e^. At the evolving horizon we have 



m{rH{t),t) 



{tV.9,,),^—-^e 



(101) 



Thus for both normalizations we have most (but not all) 
of the key properties of an outer trapping horizon at ru (t) 
when rh{t,rH{t)) > 0, i.e., when the overall mass in- 
creases in time. 



XII. DISCUSSION 

We have used the warped product formalism to inves- 
tigate the geometry of time-dependent spherically sym- 
metric spacetimes, developing relatively straightforward 
arguments for the covariant conservation of the Kodama 
vector itself and the associated Kodama flux. This con- 
struction has allowed us to construct a very general class 
of conserved fluxes appropriate to any spherically sym- 
metric spacetime. 

Furthermore, we have successfully used the Kodama 
vector field, plus the Clebsch decomposition, to obtain a 
preferred Kodama time coordinate, and have then pro- 
ceeded to construct a geometrically preferred coordinate 
system for describing spherically symmetric time depen- 
dent spacetimes. The resulting metric is one of the most 
simple forms of the metric of a spherically symmetric 
spacetime — a diagonal metric in Schwarzchild curva- 
ture coordinates. With these coordinates there are very 
simple physical interpretations for both the Hawking- 
Israel (Hernandez-Misner/ Misner-Sharp) and Brown- 
York quasi-local masses. Although the definition of sur- 
face gravity remains somewhat ambiguous, by using the 
Kodama time as an integral part of the construction we 
have identified some very good geometrically preferred 
candidates that are compatible with known results in the 
static limit. 
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